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ABSTRACT 

We report the discovery of a new “changing-look” quasar, SDSS J101152.98-1-544206.4, 
through repeat spectroscopy from the Time Domain Spectroscopic Survey. This is an 
addition to a small but growing set of quasars whose blue continua and broad optical 
emission lines have been observed to decline by a large factor on a time scale of 
approximately a decade. The 5100 A monochromatic continuum luminosity of this 
quasar drops by a factor of > 9.8 in a rest-frame time interval of < 9.7 years, while 
the broad Hct luminosity drops by a factor of 55 in the same amount of time. The 
width of the broad Ha line increases in the dim state such that the black hole mass 
derived from the appropriate single-epoch scaling relation agrees between the two 
epochs within a factor of 3. The fluxes of the narrow emission lines do not appear to 
change between epochs. The light curve obtained by the Catalina Sky Survey suggests 
that the transition occurs within a rest-frame time interval of approximately 500 days. 
We examine three possible mechanisms for this transition suggested in the recent 
literature. An abrupt change in the reddening towards the central engine is disfavored 
by the substantial difference between the timescale to obscure the central engine and 
the observed timescale of the transition. A decaying tidal disruption flare is consistent 
with the decay rate of the light curve but not with the prolonged bright state preceding 
the decay; nor can this scenario provide the power required by the luminosities of the 
emission lines. An abrupt drop in the accretion rate onto the supermassive black hole 
appears to be the most plausible explanation for the rapid dimming. 

Key words: galaxies: active - quasars: general - accretion, accretion discs. 


1 INTRODUCTION 

The dichotomy between Type 1 and Type 2 active galactic 
nuclei (AGN) is among the most studied in astrophysics (for 
a review, see Netzer 2015). The classes are empirically deter- 
* E-mail: runnoejc@psu.edu mined; Type 1 AGN are characterized by having both broad 

t Hubble Fellow (on the order of 10® km s“®) as well as narrow (on the or- 
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der of 10^ km s“^) emission lines, whereas only the narrow 
lines are present in the spectra of Type 2 AGN. The in¬ 
termediate Type 1.8 and 1.9 classifications describe objects 
with weak or absent broad H/3 emission, respectively, but a 
broad Ha line (Osterbrock 1981). In the simplest scenario, 
these objects can all be unified by orientation (Antonucci 
1993); Type 2 AGN are viewed edge-on, with an equato¬ 
rial dusty structure obscuring the line of sight down to the 
broad-line region (BLR), while Type 1 objects have an un¬ 
obscured, face-on view. Early support for this scheme came 
from ground-breaking views of broad emission lines in Type 
2 objects observed in scattered light (Antonucci & Miller 
1985). In this picture, intermediate types are explained 
by partial obscuration or reddening by optically thin dust 
(Stern & Laor 2012). In all cases, the narrow emission lines 
originate farther from the central engine in the narrow-line 
region (NLR), and are not hidden. 

While the simplicity of the above orientation-based 
scheme is compelling, many have suggested that it does 
not fully explain the Type 1-Type 2 dichotomy. An alter¬ 
native possibility is that the link between classifications 
might be evolutionary (e.g, Penston & Perez 1984). Based 
on NGG 4151 (see also Lyutyj et al. 1984) and 3G 390.3, 
which were both initially observed to have broad emission 
lines that later disappeared. Type 2 AGNs may be Type Is 
that have “turned off”. Variations along this line invoke age, 
environment, or luminosity effects and are summarized by 
Gaskell (2014) (see also, Villarroel & Korn 2014). 

Objects that transition between Type 1 and 2 pro¬ 
vide a unique opportunity to study the evolutionary 
side of unification. Although not common, the number 
of known objects in this class has been growing. The 
term “changing-look” was initially used to describe X- 
ray AGN with variable obscuration (e.g.. Matt et al. 2003; 
Puccetti et al. 2007; Bianchi et al. 2009; Risaliti et al. 2009; 
Marchese et al. 2012), but has lately widened to include ob¬ 
jects that transition between Types 1 and 2 (an optical clas¬ 
sification). There are a handful of objects like NGG 4151 and 
3G 390.3, where the broad emission lines and quasar contin¬ 
uum fade or disappear completely below the detection limit 
(Collin-Souffrin et al. 1973; Tohline & Osterbrock 1976; 
Sanmartim et al. 2014; Denney et al. 2014; Barth et al. 
2015) as well as instances where the transition goes in the 
other direction and the quasar “turns on” (Cohen et al. 
1986; Storchi-Bergmann et al. 1993; Aretxaga et al. 1999; 
Eracleous & Halpern 2001; Shappee et al. 2014). The 
timescales for the transition are typically found to be on the 
order of years, and obscuration or changes in accretion rate 
are often the preferred explanation. However, these objects 
have all been at notably low redshift and luminosity. 

LaMassa et al. (2015) report the discovery of the first 
changing-look quasar, the most luminous and distant ob¬ 
served to date, with a redshift oi z = 0.31. SDSS 

J015957.64+003310.5 (hereafter J0159-f0033), is an X-ray 
selected AGN from the Stripe 82X survey (LaMassa et al. 
2013a,b) that transitioned from a Type 1 to a Type 1.9 
within a ~ 9 year period. Initially observed in the first 
data release of the Sloan Digital Sky Survey (SDSS DRl, 
Abazajian et al. 2003), J0159+0033 appears as a typical. 
Type 1 broad-line quasar. In 2010 it was re-observed 
in the Baryon Oscillation Spectroscopic Survey (BOSS, 
Dawson et al. 2013) of the Sloan Digital Sky Survey - HI 


(SDSS-HI, Eisenstein et al. 2011), revealing a galaxy spec¬ 
trum with only a weak quasar continuum and residual broad 
Ha. Archival X-ray observations of its bright and dim states 
show that the dimming in the X-ray flux is consistent with 
the optical continuum, with no strong evidence for obscura¬ 
tion in either band. Although variable obscuration cannot 
be emphatically ruled out, LaMassa et al. (2015) favor an 
interpretation where the changes in J0159-I-0033 are due to 
a decrease in accretion power. 

The discovery of luminous analogs of the nearby 
changing-look AGNs is significant because it demonstrates 
that such abrupt transitions can occur in quasars where the 
black holes are more massive and the physical time scales 
correspondingly longer than their low-redshift counterparts. 
Motivated by the above considerations, Ruan et al. (submit¬ 
ted) performed an archival search of the SDSS DR12 spec¬ 
troscopic database, in an effort to find more such quasars. 
Their criteria target objects at all redshifts that exhibit sig¬ 
nificant changes in repeat SDSS observations, either tran¬ 
sitioning from quasar to galaxy or vice versa. In addition 
to J0159-I-0033, which is recovered, they identify two new 
changing-look quasars. In both cases, the changes in the 
spectrum are best explained by changes in the luminosity 
of a pre-existing accretion flow rather than tidal disruption 
events (TDEs), primarily because the observed narrow-line 
ratios are indicative of long-term AGN activity. Thus it ap¬ 
pears that these are examples of a previously unappreciated 
mode of quasar variability: prolonged “on-” and “off-states.” 
This type of variability may correspond to the bursting be¬ 
havior found in simulations, such as those of Novak et al. 
(2011), where the accretion rate onto the black hole is inter¬ 
mittent because of the interaction of the quasar outflow with 
the gas supplied by stellar processes (see also the related dis¬ 
cussion by Schawinski et al. 2015 and Merloni et al. 2015). 
Alternatively, the variability may be the result of thermal in¬ 
stabilities that drive the accretion disk to transition between 
different physical states (Ichimaru 1977). Such thermal in¬ 
stabilities may follow from small changes in the accretion 
rate, but the resulting large, abrupt change in the luminosity 
in this case would come from the transition of the disk state, 
e.g., to an advection-dominated flow (Narayan & Yi 1994), 
rather than a sudden cessation in accretion and draining of 
the disk. 

In this work, we report the discovery of a new changing- 
look quasar, the first found via visual inspection of the Time- 
Domain Spectroscopic Survey (TDSS, Morganson et al. 
2015). The source, SDSS J101152.98-f544206.4 (hereafter 
JlOll-l-5442), is a 2 : = 0.246 quasar initially observed 
for spectroscopy by the SDSS in early 2003 that ap¬ 
peared in the third, fifth, and seventh SDSS quasar catalogs 
(Schneider et al. 2005, 2007, 2010). It was targeted again 
by TDSS in the Few Epoch Spectroscopy (FES) program 
to study variability of the broad Balmer emission lines in a 
sample of normal, bright, low-redshift quasars and a second 
spectrum was obtained in early 2015. The main goals of this 
work are to describe the behavior of JlOll+5442 and set the 
stage for a systematic search of TDSS for a sample of chang¬ 
ing look AGN. After collecting existing data for this object 
(Section 2), we evaluate several physical scenarios that might 
explain the changes in JlOll+5442 (Section 3). We deter¬ 
mine that the changes in spectral state are driven by changes 
in the accretion power, ruling out variable obscuration and 
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TDEs. These results are summarized and discussed in the 
context of previous work in Section 4. Throughout this work, 
we adopt the same cosmology used by LaMassa et al. (2015) 
of Ho = 70 km s“^ Mpc“^, Ha = 0.73, and Q.m = 0.27, 
which leads to a luminosity distance of Dl = 1.244 Gpc 
and a distance modulus of p = 40.47 mag for JlOll+5442. 

cS 

a 

2 DATA AND MEASUREMENTS ^ 

2.1 New and archival data 

For this work, we collect the existing photometry and spec¬ 
troscopy of JlOll-l-5442. The bright-state spectrum was 
taken on January 13, 2003 through a 3V0 diameter fiber 
and covers the range 3800 — 9200 A with a resolution of 
R ~ 2000 (York et al. 2000; Smee et al. 2013) and was cali¬ 
brated with the v5_3_12 version of the spectroscopic pipeline. 
Imaging and photometry are also available from the SDSS, 
taken almost a year earlier on February 13, 2002. The SDSS 
image of JlOll+5442 shows that the host galaxy is extended ,—, 
with a visually asymmetric, potentially disturbed, morphol- ^ 
ogy. Because the source appears extended, we adopt the S 
Petrosian magnitudes, which use a large aperture (a 3”6 
diameter in the r band, corresponding to 21 kpc) to cap¬ 
ture virtually all of the light. These ugriz magnitudes are 
18.32 ± 0.01, 18.20 ± 0.01, 17.90 ± 0.08, 17.64 ± 0.09, and 
17.52 ± 0.03. The corresponding absolute magnitude in the 
i band is Mi = —22.82, thus it was included in the SDSS 
DRIO quasar catalog (Schneider et al. 2010). From the ugriz 
colors and the transformations of Jester et al. (2005) appro¬ 
priate for quasars of z < 2.1, we obtain V = 18.01 ± 0.01 
and R = 17.64 ± 0.01. The object was re-observed during 
the TDSS, a subprogram of the Extended Baryon Oscilla¬ 
tion Spectroscopic Survey (eBOSS), on February 20, 2015. 
Specifically, JlOll-l-5442 was targeted by an FES program 
to study the variability of the Balmer lines in high signal- 
to-noise ratio quasar spectra. The TDSS spectroscopic fiber 
had a diameter 2'.'0 and was centered on a position within 
fUf of the SDSS fiber position. The new TDSS spectrum 
covers the range 3600 — 10000 A at i? 2000 and was cali¬ 
brated with the v5_7_8 version of the spectroscopic pipeline. 

We apply two corrections to the spectra. To correct for dust 
extinction in the Milky Way, we deredden the spectra using 
a Seaton extinction law (Seaton 1979) with E(B — V) = 0.01 
taken from the maps of Schlegel et al. (1998). We then shift 
the spectrum to the rest frame using the SDSS pipeline red- 
shift of z = 0.246. 

In addition to the SDSS photometry in the dim state, 
there are infrared data available from the Wide-field In¬ 
frared Survey Explorer (WISE, Wright et al. 2010) ALL¬ 
WISE multi-epoch photometry catalog. JlOll-l-5442 was ob¬ 
served to have magnitudes of VFl = 13.99 ± 0.03, W2 = 
12.86 ± 0.06, W3 = 9.7 ± 0.1, and IF4 = 7.2 ± 0.2 on April 
28, 2010 and Wl = 14.19 ± 0.09 and W2 = 13.08 ± 0.08 
on November 5, 2010, where we list standard deviation of 
observations within days of each other as the magnitude un¬ 
certainty. These data do not appear to offer any additional 
information on the the variability properties of this object. 
JlOll-l-5442 has not been observed by the Chandra X-ray 
Observatory or XMM-Newton and it was not detected in the 
ROSAT all-sky survey (Voges et al. 1999). 



2002 2004 2006 2008 2010 2012 2014 
Year 


Figure 1. Top: The light curve of JlOlH-5442, including pho¬ 
tometry from the CSS (solid red circles), LINEAR (solid blue 
triangles), and the SDSS (open black square). The V magnitude 
inferred from the TDSS spectrum, observed on MJD 57073, is 
19.350 ib 0.002 and is not plotted here . This value is almost cer¬ 
tainly a lower limit on the true brightness of the AGN and host, 
as discussed in Section 2. The horizontal bar represents the likely 
dimming timescale, ~500 days in the rest frame. Bottom: The 
median light curve, where each point represents the median of 
the observations in one year bins and the error bars represent the 
rms scatter in each bin. 


31011+5442 was observed regularly between December 
2002 and June 2008 as part of the MIT Lincoln Laboratory 
Lincoln Near-Earth Asteroid Research (LINEAR) survey 
(Stokes et al. 2000). This survey is performed using a single, 
broad filter, so we adopt the transformation of Sesar et al. 
(2011) to recalibrate the LINEAR photometry: 


g — Qlinear — 0.0574 — 0.004 {g — i) 

+ 0.056 {g - if - 0.052 {g - if + 0.0262 {g - i)'‘(l) 

To apply the above equation, we use the SDSS g — i color 
and convert the light curve to the V band using the SDSS g 
magnitude and the Jester et al. (2005) transformations. The 
LINEAR light curve is shown in Figure 1 and discussed in 
more detail below. 

Overlapping the LINEAR coverage, JlOll+5442 was 
also observed between January 2006 and May 2013 as part 
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of the Catalina Sky Survey (CSS, Drake et al. 2009). Al¬ 
though the observations for this survey are made in white 
(i.e. unfiltered) light, the light curves are reported in D-band 
magnitudes. We further refine the reported magnitudes us¬ 
ing the listed transformation as a function of the V — R color: 

V = Vcss + 0.91 (V - Rf + 0.04. (2) 

This transformation was selected from the three available 
because at the redshift of this quasar neither filter includes 
the strongest emission line, Ha, while both filters include 
the second strongest emission line, H/3, in the spectral range 
where they overlap. Therefore, the V — R color suffers the 
least from emission-line contamination, compared to other 
colors we could have used. 

The total light curve, including the LINEAR and CSS 
photometry, for JlOll+5442 is shown in Figure 1. This fig¬ 
ure also includes the SDSS magnitude, which matches the 
the LINEAR magnitudes and early-time measurements in 
the CSS light curve fairly well. The late-time TDSS fiber 
magnitude is E = 19.350 ± 0.002, but this corresponds to a 
2'.'0 aperture, which excludes some of the light from the ex¬ 
tended emission from the host galaxy. Therefore, we regard 
this value as a lower limit on the brightness of the AGN and 
host and do not include it in Figure 1. It is not unusual to see 
variability in a quasar light curve; variability on the order 
of tenths of a magnitude is typical for quasars on timescales 
of a year (e.g., MacLeod et al. 2012). However, the system¬ 
atic decline of the light curve by half of a magnitude, as 
displayed in Figure 1, over a timescale shorter than several 
years is not typical behavior. 

From the light curve, we can estimate the timescale for 
the quasar flux to decay. During the first and last years of 
observation, the average brightness of the quasar appears 
constant. We estimate the transition timescale by visually 
determining where the light curve departs from the bright 
state and arrives at the dim state. The poor sampling and 
noise in the light curve make it difficult to determine exactly 
when the dimming occurs, but our best estimate is ^500 
days in the quasar rest frame. The relevant portion of the 
light curve is marked in the figure by a horizontal bar. 


2.2 Spectral decomposition 

We decompose the spectra interactively using the IRAF task 
SPECFIT (Kriss 1994), which fits multi-component models to 
spectra by minimizing the statistic. Components for the 
fit include a power law of the form f\ = A\~°‘ to char¬ 
acterize the quasar continuum, Gaussians to characterize 
the emission lines, optical (Boroson & Green 1992) and UV 
(Vestergaard & Wilkes 2001) Fell templates that have been 
combined into a single template that covers our full wave¬ 
length range, and Chariot & Bruzual starlight templates (see 
Gales et al. 2013). The spectral decomposition is performed 
over the entire available wavelength range (3100-8000 A in 
the dim state and 3100-7500 A in the bright state) and all 
components are fitted simultaneously. 

First, we decompose the dim-state spectrum because 
this spectrum is simpler than the bright-state spectrum. 
The starlight and narrow-line properties inferred from this 
spectrum aid us in the decomposition of the bright-state 


spectrum, as we explain below. We impose the following 
constraints on Gaussians representing narrow lines: within 
a doublet, the two emission lines are constrained to have 
identical velocity widths, a flux ratio of 3 when appropri¬ 
ate, and a constant wavelength ratio. We further group the 
lines by ionization and critical density and constrain the 
velocity widths to be the same across a group; Ha and H/1 
form the first group, [O ii] AA3726, 3729, [N ii] AA6548, 6583, 
and [Sii] AA6716, 6731 the second, and [Neiii] AA3869, 3967 
and [Olll] AA4959, 5007 the last. Finally, we also tie the 
wavelength of the narrow H/3 to that of the narrow Ha, be¬ 
cause the Ha line is clear in the spectrum and the H/I line 
from the quasar is weak enough that it is severely contam¬ 
inated by the stellar lines. In order to achieve a satisfac¬ 
tory fit, this spectrum requires two stellar populations: one 
where the age is fixed at 10 Gyr and the normalization is a 
free parameter, plus a second where both the age and nor¬ 
malization are allowed to vary. This approach is similar to 
Ganalizo & Stockton (2013), who find significant contribu¬ 
tions of intermediate plus old stellar populations in quasar 
host galaxies. The broad lines (including the Fell) in this 
spectrum are very weak; only a symmetric broad Ha profile 
is visible and it is well characterized by a single Gaussian. 
The resulting fits are shown in Figures 2 and 3. 

In the bright state, we fit the spectrum with all the com¬ 
ponents used in the dim state, plus the Fe ll template and a 
suite of Gaussians to characterize the broad lines. The broad 
Ha and H/3 lines are represented by two Gaussians each in 
order to describe their asymmetries, while all other broad 
lines are represented by single Gaussians. Since the size of 
the NLR {Rnlr ~ 2.5 kpc based on the dim-state [O ill] lu¬ 
minosity and the correlation of Hainline et al. 2013) should 
fall within the TDSS fiber, we fix all narrow-line properties 
to the best-fitting values from the dim-state decomposition. 
For the starlight templates, we fix the ages to the best-fitting 
values from the dim state but allow the normalizations to 
vary. Although we do not expect intrinsic variability from 
the host galaxy light, the size of the fibers used to obtain the 
observations were 3'.'0 and 2'.'0 in the bright and dim states, 
respectively, which creates the potential for differences be¬ 
tween the spectra. Visual inspection of the SDSS image of 
this object reveals that the host galaxy extends beyond the 
diameter of the fibers, suggesting that the larger fiber may 
include more host-galaxy light. Thus, we do not require that 
the normalization of the stellar templates be identical in the 
two states. 

The uncertainties in all parameters resulting from the 
spectral decomposition are determined via a full Monte- 
Garlo simulation of the fitting procedure. We resample each 
spectrum by selecting for each pixel a new value for the 
flux density, drawn from a Gaussian distribution that has a 
mean equal to the measured flux density in that pixel, and 
a standard deviation equal to the measured uncertainty in 
the flux density. The best-fitting parameters for the observed 
spectrum are used as an initial guess, and the synthetic spec¬ 
trum is refitted. This process is repeated 10® times each for 
the bright- and dim-state spectra to generate distributions 
of all spectral fit parameters. The standard deviations of 
these distributions are taken to be the Icr uncertainties in 
the measured properties. 

The complete model fits to the bright- and dim-state 
spectra are shown in Figure 2. In Figure 3 we show an ex- 
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panded view of three spectral windows of particular interest. 
In the latter figure we also illustrate the individual model 
components. The bright and dim state emission-line mea¬ 
surements are listed in Table 1; the continuum properties we 
describe here. In the bright state, the best-fitting power law 
has a normalization of (1410±60) x 10“^^ erg s“^ cm“^ 
at 1000 A and an exponent of 3.42 ± 0.06. In the dim state, 
the best-fitting normalization and exponent are (11 ± 2) x 
10“^^ erg s“^ cm“^ A“^ and 1.9 ± 0.3, respectively. The 
full width at half maximum (FWHM) of the lines in the 
Fell complex, which is only present in the bright state, is 
1500 ± 200 km s“^. The age of the young stellar template is 
225 ± 4 Myrs in both states. The flux ratio of the young to 
old stellar templates at 5500 A is 0.21 ± 0.01 in the bright 
state and 0.76 ± 0.01 in the dim state. We attribute the dif¬ 
ference in this ratio to the difference in seeing conditions 
and fiber diameters and interpret it to mean that the young 
stellar population is more centrally concentrated than the 
old stellar population. In the dim state, the quasar contin¬ 
uum is detected above the la level of the noise in the data 
at wavelengths shorter than approximately 6800 A and the 
broad H/3 line is not detected at all. We can place a limit 
on the dim-state broad H/3 line by asking what broad 
fluxes could be hidden in the noise in the optical part of the 
spectrum; using this criterion the broad H/3 integrated flux 
must be less than 13.7 x 10“^^ erg s“^ cm“^ at the Icr level. 

Motivated by the apparently weak power-law contin¬ 
uum in the dim state, we also test our model assumptions 
for the dim-state fit by asking whether we can obtain a sim¬ 
ilarly good decomposition without a power-law component. 
We perform a second fit to the dim-state spectrum, identi¬ 
cal to our adopted dim-state decomposition described above 
but with the power-law component removed. Visually, the 
result describes the data well. We compare the Bayesian In¬ 
formation Criterion (BIC Schwartz 1978), similar to the chi- 
squared statistic but appropriate for a comparison of models 
with a different number of free parameters, for this new ht 
and our adopted dim-state fit and find that statistically they 
describe the data equally well. The emission line components 
of the ht are essentially unchanged compared to our adopted 
ht, so the result of removing the power-law continuum is that 
the normalizations of the stellar continuum components in¬ 
crease and the age of the young stellar population decreases 
to 100 Myrs. Thus, the best-htting power-law continuum 
in the dim state and the properties derived from it (e.g., 
luminosities) should be considered upper limits. 

2.3 Inferred quasar properties 

Using the decomposed bright-state spectrum, we measure 
fundamental properties of the central black hole. In the 
bright state, we apply the 5100 A bolometric correction with 
a nonzero intercept from table 1 of Runnoe et al. (2012) 
to the monochromatic optical luminosity A La (5100 A) = 
(6.3±0.1) X 10^® erg s~^ to calculate a bolometric luminosity 
of Uoi = (6.8±0.1) X lO^"* erg s'A The broad H/3 FWHM of 
2350 ±50 km s“^, combined with the single-epoch mass scal¬ 
ing relationship from (Vestergaard & Peterson 2006), yields 
Mbh = (3.6 ± 0.2) X 10^ Mq. These estimates produce an 
Eddington ratio of L/Lstid = 0.15. We can also calculate 
the black hole mass from Ha by adopting the scaling re¬ 
lationship from (Greene et al. 2010). This calibration dif¬ 


fers slightly from that of Vestergaard & Peterson (2006), 
as discussed by Gales et al. (2013), but the difference is 
small compared to the intrinsic scatter. This approach yields 
Mbh ~ (4.7T0.1) x 10^ Mq and L/LEdd = 0.12. The uncer¬ 
tainties quoted on these mass measurements do not include 
the 0.4 — 0.5 dex scatter associated with the single-epoch 
mass scaling relationships. This additional source of error 
translates to a factor of 2.5 — 3 uncertainty in mass and 
is more representative of the uncertainty than the formal 
errors. Thus, the Ha and H/l masses are consistent with 
each other given the large uncertainties associated with the 
mass scaling relationships. Using the same relationships for 
the dim state, ALa(5100 A) = (6.3 ± 0.1) x lO'^^ erg s“^, 
Lboi = (8.4T0.1) X lO'*® erg s"\ Mbh = (1.8T0.3) x 10^ Mq, 
and L/LEdd = 0.04 (the mass measurement is for Ha as H/1 
is not detected). 

In principle, because the black hole mass is constant, the 
continuum luminosity and broad-line widths should change 
in concert such that the single-epoch mass measurements are 
consistent between the bright and dim states. LaMassa et al. 
(2015) found this to be the case for J0159±0033. Our esti¬ 
mates for the Ha masses are preserved within a factor of 
three, although not within the formal Icr uncertainties. 

To constrain the excitation mechanism of the narrow¬ 
line gas, we plot the observed [Oill]/H/3 and [Nil]/Ha 
narrow-line ratios in a diagnostic diagram (also known as a 
“BPT” diagram after Baldwin, Phillips & Terlevich 1981). 
The location of the measured ratios in this diagram can 
reveal the relative contributions to the ionizing continuum 
from the AGN and star formation. The line ratios from the 
dim-state spectrum, where the narrow-line decompositions 
as performed, are shown in Figure 4. We include the classifi¬ 
cation schemes of Kewley et al. (2001) and Kauffmann et al. 
(2003) to separate line ratios characteristic of AGN, star¬ 
forming, and composite galaxies and show as contours the 
distribution of SDSS galaxies from Tremonti et al. (2004) 
for comparison. The conclusion for J1011±5442 is clear: a 
hard AGN continuum is required to produce the observed 
narrow lines. This does not change if we adopt the alterna¬ 
tive line ratios of ]Oi] A6300/Ha or ]Sii] AA6716, 6731/Ha 
in place of ]N ll]/Ha. 


3 ANALYSIS 

In this section, we perform a series of tests on the spectra and 
LINEAR±CSS light curve in order to constrain the physi¬ 
cal mechanism responsible for the drastic changes observed 
in J1011±5442. We consider three physical scenarios for the 
observed dimming of the quasar central engine: abrupt ob¬ 
scuration of the central engine, intrinsic dimming due to a 
drop in the accretion rate (also examined as explanations for 
the dimming of J0159±0033 by LaMassa et al. 2015), and 
decay of a flare produced by the tidal disruption of a star 
by the central BH (proposed in the context of J0159±0033 
by Merloni et al. 2015). 

3.1 Obscuration of the central engine 

Following LaMassa et al. (2015), we examine whether the 
time required for a foreground object in a bound orbit 
around the BH to move in front of the continuum source 
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Rest Wavelength [A] 

Figure 2. The spectral fits for the bright-state (MJD 52652) and dim-state (MJD 57073) spectra. The data are shown in black and 
the best-fitting total model is shown in red. See Section 2 for a description of the spectral decomposition and Figure 3 to see the flux 
uncertainties. 


Table 1. Properties of strong emission lines 

Bright State 

A FWHM 

Line*^ [A] [km s“^] 

Integrated 

Flux*’ 

A 

[A] 

Dim State 
FWHM 
[km s“^] 

Integrated 

Flux*’ 

[O ll] n 

e 



3729.0±0.2 

499i7-f 

20il 

[Nelll] n 




3970.0±1.0 

1310i60f 

lOil 

H(5 b 

4110.0±1.0 

4170±60 

99±6 




Hy b 

4346.0±0.7 

2800±40 

290±3 




H/3 bl'^ 

4867.0±0.4 

1960±60 

530±40 




b2'= 

4860.0±0.9 

3800±200 

450±40 




bd 

4864.0±0.3 

2350±50 

980±10 



il3.7 

n 




4861.9±0.1-f 

361i7f 

nil 

[Olll] n 




5001.0±0.4 

1310i60f 

82i3 

Hel b 

5895.0±1.0 

3200±300 

200±10 




[Ol] n 




6363.0±0.5 

670i30 

5il 

[N ll] n 




6584.0±0.1 

499i7f 

50i2 

Ha bl'^ 

6568.0±0.1 

2260±20 

3430±40 

6563.5±0.1 

2720i70 

87i3 

b2'= 

6580.0±1.0 

10200±400 

1340±30 




bd 

6571.0±0.0 

2420±30 

4770±30 

6563.5±0.1 

2720i70 

87i3 

n 




6564.0i0.l-f 

361i7f 

34il 

[S ll] nl 




6717.0i0.6-f 

499i7f 

23il 

n2 




6732.0i0.6-f 

499i7f 

17il 


^ b = broad component (number identifies Gaussian components summed to produce total 
line profile), n = narrow component. Wavelengths for lines that are in a doublet are for the 
long-wavelength line. 

^ Integrated fluxes are in units of 10“^^ erg s“^ cm“^. 

^ One of two Gaussian components. 

^ Total profile, the sum of two Gaussian components. 

^ The bright-state narrow lines parameters are fixed to the dim-state values. 

^ This value is constrained during the spectral decomposition. See Section 2.2 for a descrip¬ 
tion of the constraints. 


and BLR can be short enough to match the observed tran¬ 
sition time. We adopt equation 4 of LaMassa et al. (2015), 


= 0.07 


^ orb 


1 It — day 


3/2 




- 1 / 2 ^ 


^ orb 


yr, (3) 


which gives the time for an object at an orbital radius 
Verb to traverse a region extending to radius Vsre around 


a black hole of mass Mg — M/10 ^Mq. The most conser¬ 
vative case, which minimizes the crossing time, corresponds 
to the minimum plausible value of Vsre and a value of Vorb 
that is determined graphically. Thus, we begin by assuming 
that Vsre is the radius derived from the “Clean2-|-ExtCorr” 
radius-luminosity relation {Rblr, Bentz et al. 2013), i.e., 
an emissivity-weighted radius that encompasses the bulk of 
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Figure 3. The spectral decomposition in three wavelength windows including Ha (top), H/3 (middle), and strong stellar absorption lines 
(bottom). The data (black), uncertainties (gray), best-fitting model (red), old and young stellar templates (yellow, orange), power-law 
continuum (green), and emission lines (blue, see Table 1 for a list) are displayed for the bright state (MJD 52652) on the left and the 
dim state (MJD 57073) on the right. 
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log([N II] X6583 / Ha X6563) 


Figure 4. Diagnostic emission line ratios for JlOll+5442. The 
narrow lines in the bright-state spectral decomposition were 
fixed to the dim-state values, so only one position is obtained 
on the BPT diagram. The BPT identification schemes from 
Kewley et al. (2001) and Kauffmann et al. (2003) are shown as 
dotted and dashed lines, respectively. For comparison, gray con¬ 
tours show the properties of all SDSS galaxies in Tremonti et al. 
(2004). The emission-line ratios in JlOll-l-5442 are consistent 
with a pure AGN ionizing continuum. 


the broad-line emission. For JlOll-l-5442 Rblr = 28 light- 
days; assuming that Vsrc = Rblr, we estimate a minimum 
crossing time of tcross ~ 21 years with Vorb = 34 light days. 
This minimizes the crossing time since we have effectively 
assumed that the obscuring object is orbiting within the 
BLR and is large enough to obscure its brightest parts. The 
absolute lower limit on the crossing time is toross ~ 9 years, 
calculated by subtracting the 0.13 dex scatter in the R — L 
relationship from log{RBLR), increasing the mass by a fac¬ 
tor of 2.5 to account for the uncertainty in the mass scal¬ 
ing relationship, and recalculating the crossing time un¬ 
der the same assumptions about Vsrc and Vorb- However, 
more realistically, the obscuring object should orbit out¬ 
side the BLR, whose outer radius is several times larger 
than Rblr- Reverberation mapping of infrared emission 
from illuminated dust indicates that the BLR is contained 
within ~ 3 Rblr (e.g., Barth et al. 2013; Suganuma et al. 
2006; Baldwin et al. 2004, see LaMassa et al. 2015 for a 
more detailed justification of this assumption), so we set 
Lorb = rare = 3Rblr- TMs leads to a crossing time of 
toross ~ 31 years. To obscure all of the BLR, or for screens 
located farther from the source, crossing times can easily 
reach hundreds of years. In all cases, these timescales are 
far too long to explain the dimming seen here. 

If the power-law continuum is detected, we can carry 
out a complementary test by comparing the ratio of fluxes 
from the dim- and bright-state spectra at specific wave¬ 
lengths and asking whether they are all consistent with the 
same amount of reddening. The ratio of the Ha flux between 
the dim and bright states is 0.018 ± 0.001, which requires 
E{B — V) = 1.744 in order to be explained by extinction. 
The flux ratio of Ha/H/3 changes by > 1.3 from the bright to 


dim state, implying E(B — V)> 0.248, which is consistent 
with the decline in the Ha flux. The ratio of the quasar con- 
tinua from the two epochs at 6563 A is 0.15 ± 0.08. which is 
not consistent with the amount of dimming measured at the 
same wavelengths in the Ha line. The H)I line and nearby 
continuum should also demonstrate the same amount of at¬ 
tenuation, but this line ratio is < 0.014 whereas the 5100 A 
continuum is 0.100 ± 0.003. Finally, at 3500 A the amount 
of reddening needed to attenuate Ha produces a continuum 
flux that is about a factor of 100 too weak. Thus, if the 
continuum is detected from the quasar the relative dimming 
of the emission lines and the continuum is not consistent 
with reddening. If, on the other hand, the power-law con¬ 
tinuum is an upper limit, this discrepancy is alleviated and 
the real problem for the reddening scenario is the timescale 
for obscuration to set in. 

Additional evidence against the reddening scenario can 
be found in the broad Ha line width; the line profile changes 
such that Ha is broader in the dim state, suggesting that 
it is emitted nearer to the central engine which would pre¬ 
sumably suffer from the most extinction. Based on all these 
calculations, we disfavor reddening by a foreground screen 
as an explanation of the changes in this object. 

3.2 Intrinsic dimming of the continuum source 

An alternative explanation for the observed variability is 
that a rapid change in accretion rate is responsible for 
the dimming in the continuum and broad emission lines. 
We evaluate whether this might be possible by calculat¬ 
ing the accretion timescale for the inner part of the accre¬ 
tion disk where most of the ionizing radiation is generated 
{R < 10 Rg, where Rg = GMbh/R is the gravitational ra¬ 
dius). We follow equation 5 of LaMassa et al. (2015) for the 
inflow timescale (i.e., the time for a parcel of gas to move 
radially from a given radius in the accretion disk to the 
center, also known as the viscous or radial drift timescale), 
adopting their fiducial value for the viscosity parameter 
(Shakura & Sunyaev 1973) but substituting our measured 
values of Ms = 0.36 and \Edd = 0.15. This gives an inflow 
time of tinfi = 0.25 years, or equivalently three months. This 
timescale is short enough that a drop in accretion can eas¬ 
ily explain the changes in the continuum that are observed 
in JlOll-l-5442. The inflow time scale we estimate here is 
considerably shorter than what LaMassa et al. (2015) esti¬ 
mated for J0159-I-0033 because the two quasars have con¬ 
siderably different masses and accretion rates (in the bright 
state, JlOll-l-5442 has Ms = 0.36 and Xsdd = 0.15 versus 
Mg = 1.7 and \Edd = 0.04 for J1059-b0033). 

It is interesting to note that the continuum throughout 
the optical band dimmed within approximately 500 days, 
as well as the far-UV continuum at ~ 1 Ry, which is re¬ 
sponsible for the ionization of Hydrogen and the produc¬ 
tion of the Balmer lines. As noted by LaMassa et al. (2015), 
these changes may be driven by a drop in the luminosity of 
the inner accretion disk even though the optical continuum 
is produced in a region of the accretion disk that is much 
larger than 10 Rg. This because the optical continuum may 
result from illumination from the inner disk, as suggested 
by the results of reverberation studies (e.g., Krolik et al. 
1991; Cackett et al. 2007). Similarly, since the “radius of 
the broad-line region,” Rblr, is only 28 light-days, the lu- 
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minosity of the broad Balmer emission lines will track the 
decline of the Lyman continuum with a lag that is less than 
a month. Therefore, the transformation of JlOl 1+5442 can 
be accomplished by a rapid drop in the far-UV continuum 
luminosity resulting from an abrupt decline in the accretion 
rate. The mechanism behind this change in the accretion 
rate, however, is unknown. Tantalizing as this scenario may 
be, it must remain tentative until a plausible mechanism for 
abrupt changes in accretion rate is identified. 

3.3 Tidal disruption events 

A decaying flare from a tidal disruption event (TDE) was 
proposed as a possible mechanism for the changes ob¬ 
served in J0159+0033 by Merloni et al. (2015). Specifically, 
J1059+0033 was initially observed in a lowly accreting state, 
but brightened dramatically before spectroscopy revealed 
the presence of broad lines, which later disappeared. We 
use the LINEAR+CSS light curve to evaluate if the TDE 
scenario can explain JlOl 1+5442. In this case, there is no 
dramatic brightening of the light curve, so the entire event 
is not covered. Rather, the hypothetical TDE is already un¬ 
derway at the start of our observations so we resolve only 
part of the peak of the light curve followed by its decline. 

For a more quantitative investigation, our first step is to 
derive the light curve of the variable component by subtract¬ 
ing the flux of the host galaxy. Although the TDSS spectrum 
provides a good measure of the host emission, the 2”0 aper¬ 
ture likely misses some of the flux. Therefore, instead of us¬ 
ing the TDSS spectrum, we assume that the faintest point in 
the CSS light curve provides a limit on the host flux, which 
we subtract from all points in the light curve. 

Analytic calculations (Phinney 1989) and numerical 
modeling (e.g., Evans & Kochanek 1989; Ayal et al. 2000) 
of TDEs predict that the light curve of a TDE will decline as 
where a = —5/3. In principle, this decay law that 
is traditionally associated with TDEs actually describes the 
rate at which material will return to the vicinity of the black 
hole, which is likely not equal to the accretion rate. In fact, 
Lodato & Rossi (2011) predict different power-law indices 
for light curves in different electromagnetic bands. Nonethe¬ 
less, the majority of the tidal disruption events claimed to 
date do exhibit this decay law in their light curves (e.g., 
Gezari et al. 2009). To test whether the declining light curve 
of JlOll+5442 is consistent with this expectation, we write 
the evolution of its E-band magnitude as 

V{t) = Vo + alog(^^^y (4) 

where to indicates the beginning of the decline in the light 
curve, observationally constrained to be between 2009.4 and 
2009.9. To fit the light curve, we adopt to = 2009.65 and use 
the range in to to determine the uncertainty on a. We find 
a = —1.31q 7, which is formally consistent with —5/3; both 
curves are presented with the host-subtracted light curve 
in Figure 5. We note that, for some specific values of to 
within 2009.4-2009.9, it is possible to obtain slopes that are 
inconsistent with —5/3. 

We can also calculate the amount of matter accreted 
by the black hole during the bright state and compare it to 
reasonable values for stars that might be tidally disrupted. 
That is to say, the amount of matter accreted during the 



Year 


Figure 5. The Catalina Sky Survey light curve of the variable 
emission with the host subtracted. Red points indicate data that 
were included in fitting the TDE model, while other points are 
shown with gray error bars only. The gray region indicates the 
observational constraints on to and the dashed gray line is the 
adopted value. The solid black line represents the best-fitting 
TDE model, where a = —I.S/Iq®, and the dotted line shows 
the same model with the slope fixed at a = —5/3. 


high state cannot exceed the mass of the disrupted star. As¬ 
suming an efficiency of 77 = 0.1, and adopting the measured 
bolometric luminosity from the high state, the black hole 
accreted 0.7 Mq in the ~6 years that it was in the bright 
state. 

Although the above tests appear to support the TDE 
interpretation, the overall shape of the light curve does not. 
The most recent simulations for feeding black holes with 
tidal disruption of stars show that the light curve of a TDE 
has a very narrow peak Guillochon & Ramirez-Ruiz (2013). 
From their figure 5, a timescale of a few months is a generous 
upper limit on the bright peak in the light curve of a TDE. 
Similarly, the observed light curves of strong TDE candi¬ 
dates (e.g., Gezari et al. 2009) remain in the bright state 
only for a very short time, timescales on the order of months. 
In contrast, JlOll+5442 is in the bright state for at at least 
years in the rest frame. 

The bright-state emission lines are difficult to explain 
in the TDE scenario. The observed broad emission lines sug¬ 
gest a BLR that is too extended and massive to have been 
deposited by the recent tidal disruption of a single star. The 
narrow lines, which do not change, also present challenges, 
requiring a hard ionizing continuum and gas at even larger 
radii. Merloni et al. (2015) discuss these points and relevant 
calculations in detail, but the issue remains that, while a 
TDE can illuminate an existing BLR, it cannot distribute 
the necessary gas around the black hole or produce narrow 
lines coincident with the peak of the light curve. Based on all 
these lines of evidence, we therefore reject the TDE scenario 
to describe JlOll+5442. 
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4 SUMMARY AND DISCUSSION 

In this work, we report the discovery and analysis of the a = 
0.246 changing-look quasar, SDSS J101152.98-1-544206.4. In 
the bright state, observed in early 2003 as part of the 
SDSS, JlOll-l-5442 Is listed in the SDSS quasar catalogs up 
to DR7 and the optical spectrum of JlOll-l-5442 exhibits 
the broad emission lines and a blue continuum of a typ¬ 
ical broad-lined quasar. From this spectrum, we determine 
Lioi = (6.8±0.1) X10'^'* erg s'S Mbh = (3.6±0.2) x lO'^ Mq, 
and \Edd = 0.15. The TDSS targeted the object for obser¬ 
vations in early 2015 and found It in a dim state: the quasar 
continuum was undetected at the la level at wavelengths 
longer than ~ 6800 A, and of the broad emission lines, only 
Ha was still present, but at a much weaker level than was 
observed during the bright state in 2003. The 2006—2013 
CSS light curve reveals that the transition likely occurred 
over ~ 500 days in the rest frame around 2010. 

We test several physical scenarios to explain the dra¬ 
matic changes in J1011-|-5442’s spectrum. In some cases 
(typically in Seyfert 1.8-1.9 galaxies), variable obscuration 
of the central engine can explain changing-look AGN, but 
this mechanism is disfavored in this case. The timescale for 
an intervening cloud of dust to obscure the broad-line re¬ 
gion is too long, on the order of 30 years, compared to the 
measured timescale of ~ 500 days. An additional test of 
the obscuration scenario would be to ask whether the ob¬ 
served change in the quasar continuum is consistent with the 
amount of reddening needed to produce the change in Ha. 
We attempt this test, but the results are inconclusive be¬ 
cause the quasar continuum is not robustly detected in the 
optical. Hubble Space Telescope observations in the ultra¬ 
violet, where the power-law continuum rises, would be re¬ 
quired to perform this test conclusively. The observations 
of JlOll-l-5442 cannot be explained by a TDE. Although 
the decay rate of the light curve is consistent with theoret¬ 
ical expectations for TDEs and the 0.7 Mq accreted while 
JlOll-l-5442 was in the bright state could have been de¬ 
posited by disruption of a star, the bright state has a mini¬ 
mum duration that is far longer than any putative observed 
TDE or theoretical TDE model. An abrupt change in the 
accretion rate can drive the changing look of JlOll-l-5442. 
The timescale for material to be transported through the 
inner accretion disk and for the disk to respond to changes 
in the accretion rate is on the order of months, consistent 
with the observed change, making this a viable process. 

In the context of recently discovered changing-look 
quasars (LaMassa et al. 2015, Ruan et ah, submitted; Gales 
et al. in preparation), JlOll-l-5442 has similar luminosity, 
redshift, and transition timescale. All recently-discovered 
changing-look quasars cluster in the upper right corner of 
the luminosity-redshift diagram of LaMassa et al. (2015, see 
their Figure 1). JlOll-l-5442 does, however, have a notably 
smaller black hole mass and larger Eddington ratio com¬ 
pared to J0159-I-0033, which lead to a much shorter inflow 
time scale in the inner accretion disk (months compared 
to decades for J0159-I-0033). As a result, a large, abrupt 
fluctuation in the accretion rate is a viable explanation of 
the variability of JlOll-l-5442, whereas with J0159-I-0033 
there is some potential for ambiguity. Elitzur et al. (2014) 
describe an accretion-rate driven paradigm for the evolu¬ 
tion between spectral types. This is based on the sample of 


intermediate-type AGNs from Stern & Laor (2012), where 
a decreasing fraction of the bolometric luminosity is con¬ 
verted to broad-line emission as the accretion rate, charac¬ 
terized by Lhoi/M^^, drops. J1059-I-0033 is consistent with 
the expected values of Lboi/Mg^^ for this sequence but it 
is not possible to ascribe the spectral changes to changes 
in the accretion rate with complete certainty because the 
inflow timescale is so long, whereas JlOll-l-5442 is a fairly 
unambiguous case because of the shorter inflow timescale. 
Thus, the changing look-quasars appear to be more lumi¬ 
nous analogs of some of the Seyfert 1.8-1.9 galaxies stud¬ 
ied by Goodrich (1989, 1990, 1995) and later Trippe et al. 
(2010). As Goodrich and Trippe et al. conclude from a va¬ 
riety of observational tests, the variability of the contin¬ 
uum and broad Balmer lines of some objects could be at¬ 
tributed to intrinsic variations while in others it was caused 
by changes in the reddening. 

Our understanding of changing-look quasars will bene¬ 
fit substantially from continued optical spectroscopic mon¬ 
itoring as well as X-ray observations. X-ray spectroscopy 
was important in rejecting the obscuration scenario for 
J0159-I-0033 and can serve the same purpose for the other 
changing-look quasars. Spectroscopic monitoring can serve 
at least two purposes. First, it can reveal whether the quasar 
central engine turns on again because X-ray emission orig¬ 
inates near the black hole and acts as an instantaneous 
tracer of accretion. As we discover and monitor larger num¬ 
bers of such objects we will be able to constrain the dura¬ 
tion of the “off-state.” Moreover, if the quasar returns, this 
would be a strong indication that the previous “on-state” 
was not caused by a TDE since the rate of such events is ex¬ 
pected to be ~ 10“®-10“^ galaxy“^ (both on theoret¬ 

ical and observational grounds, e.g., Merritt & Poon 2004; 
Donley et al. 2002; Gezari et al. 2008; Maksym et al. 2010). 
Second, as LaMassa et al. (2015) remarked, we may observe 
a decline in the narrow emission lines, starting potentially as 
early as several years after the decline of the continuum. The 
precise onset of the response will depend on the structure 
of the NLR in JlOll+5442, but can happen more quickly 
than imagined given that the traditional kiloparsec-scale ra¬ 
dius of the NLR describes the outer extent rather than the 
innermost radius of the gas (see discussions on narrow-line 
variability in e.g., LaMassa et al. 2015; Denney et al. 2014). 
This would afford an additional test of the intrinsic dimming 
hypothesis and it would allow us to probe the structure of 
the narrow-line region. There are, in fact, specific predic¬ 
tions of how the narrow-lines would respond to a sharp drop 
in the ionizing continuum, based on time-dependent pho¬ 
toionization calculations by Binette & Robinson (1987) and 
Eracleous et al. (1995). Among the noteworthy predictions 
are the rapid response of the optical [O ill] doublet (the line 
flux tracks the decline of the ionizing continuum with a lag 
of only a few years) and the very slow response of the Balmer 
lines (they remain fairly steady for decades to centuries after 
the ionizing continuum has faded). 

The TDSS offers a promising way of discovering sub¬ 
stantial numbers of changing-look quasars because it will 
revisit several thousand objects with previous spectra from 
the SDSS. The two changing look quasars discovered within 
the first year of the TDSS may be the tip of the iceberg. 
By design, however, the TDSS is unlikely to find transitions 
from the “off” to the “on” state, which are just as inter- 
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esting and valuable in their own right, because the TDSS 
targets are known AGN. 
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